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1. Introduction 

Human biomagnetism has been an important area of fundamental research and 
medicine since the first low noise studies by Cohen et al. (1970). Magnetic signals 
can provide complimentary or unique information in several applications, such as 
magnetoencephalography (MEG) (Hamalainen et al. 1993) and magnetocardiography 
(MCG) (Wakai et al. 2000). For example, high quality foetal MCG (fMCG) signals can 
direct in utero diagnosis and therapy of cardiac arrhythmia such as supraventricular 
tachycardia and torsades de pointes (Wakai et al. 2003, Cuneo et al. 2003). 

Typically, these studies use superconducting quantum interference devices 
(SQUIDs) for magnetic field detection, usually operated in arrays of matched 
gradiometers with 30 or more channels. SQUIDs have achieved a noise limit of around 
lfTHz -1 / 2 (Robbes 2006), but they are expensive and require liquid 4 He for the low 
T c systems for the high signal-to-noise ratio (SNR) needed for demanding biomagnetic 
applications like fMCG (Li et al. 2004). 

Atomic magnetometers have recently shown promise for use in biomagnetism. 
Atomic vapor magnetometers operating in the spin exchange relaxation free (SERF) 
regime (Allred et al. 2002) have recently surpassed the sensitivity of SQUID systems, 
achieving 160aTHz _1//2 sensitivities (Dang et al. 2010). Additionally, the application of 
microfabrication techniques to the manufacture of mm-scale alkali-metal vapor cells and 
integrated fibre-coupled optics (Shah et al. 2007) allows the possibility of large arrays of 
SERF magnetometers, similar to current commercial SQUID systems. Besides offering 
increased sensitivity, atomic magnetometers are operated at temperatures between 25- 
180 C (Budker & Romalis 2007), replacing expensive cryogenics and maintenance of 
SQUID systems with simple resistive electrical heating and passive thermal insulation. 

In this paper, we report a four-channel portable atomic magnetometer array suited 
for biomagnetic applications. Each magnetometer in the array features a 6-11 fTHzT 1 / 2 
noise floor from 10-100 Hz and has adjustable channel spacing and orientation. All 
channels are self-contained, including local tri-axial nulling coils, in principle allowing 
the channels to be placed in an arbitrarily oriented non-planar geometry. Using this 
array in a magnetically shielded room (MSR), we have made high quality adult MCG 
measurements to demonstrate the efficacy of our device for use with human subjects. 

Previously, adult MEG has been measured using a similar sized fibre coupled atomic 
magnetometer operating in a different detection mode (Johnson et al. 2010). Single 
channel sensitivity was comparable to ours, but quadrant photodiode detection was 
used to construct gradiometer signals, rather than the use of multiple magnetometers. 
A similar technique using a photodiode array was used previously (Xia et al. 2006), 
but with a large non-portable setup. In both cases, the separation between photodiode 
elements is small compared the the signal source depth of an MCG measurement. In 
this regime, the individual sensors will measure approximately the same signal, and will 
not operate as truly multichannel devices on their own. 

Adult MCG studies have also been performed. Bison et al. (2009) used an array 
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of 25 room temperature Cs (non-SERF mode) magnetometers as a 19-channel set of 
2nd order gradiometers to map adult MCG in an eddy-current shield. They were able 
to suppress the noise by a factor of 1000 using this technique, but were limited to 
300 fT Hz -1 / 2 in this configuration. 

Recently, Knappe et al. (2010) used a single channel microfabricated mm-scale 
alkali-metal vapor cell to measure adult MCG in an MSR. Their signals compare well 
with SQUIDs simultaneously measuring the MCG. The measurements were performed 
with a sensitivity 100-200 fTHz" 1 / 2 but in a much quieter shielded room (Bork 
et al. 2000). Results from a single magnetometer were presented. 

The outline of the rest of the paper is as follows: in section [2j we summarize 
the theory of operation of our SERF magnetometer. Section [3] presents our apparatus 
and design features. Section H] details the operation and characterization procedure of 
our array. In section [5j we show a sample adult MCG, compare our measured QRS 
amplitudes from 13 subjects with previous results in the literature, and demonstrate 
the acquisition of much smaller test signals applied with a phantom. Section [6] provides 
an analysis of the noise sources limiting our sensitivity. Finally, section [7] concludes with 
directions for future work towards the goal of obtaining fMCG signals. 

2. Theory 

The magnetometer presented here operates similarly to previous SERF magnetometers 
such as Alfred et al. (2002) and Li et al. (2006). Optical pumping creates a nonzero 
average electron spin along the propagation direction of the pumping laser. Orthogonal 
magnetic fields cause the population spin to precess at the slowed down Larmor 
frequency (Happer Sz Tang 1973). Spin- relaxation collisions and the absorption of laser 
photons interrupt the precession. The balance between optical pumping, spin relaxation, 
and Larmor precession determines the dynamics of the atomic polarisation and the 
performance of the magnetometer. A probe laser detects the orthogonal components of 
the electron spin. 

Figure []Ja) and (b) show the optical setup and 87 Rb energy level diagram, 
respectively. A circularly polarised pump beam resonant with the "Dl" transition 
optically pumps the atoms. An orthogonal probe beam detuned from the "D2" 
transition detects atomic spin precession due to external fields using the Faraday effect 
and a balanced polarimeter. The polarisers at the output of the optical fibre convert 
laser polarisation noise into intensity noise, to which our configuration is much less 
sensitive. 

2.1. Polarisation equation of motion and steady state solution 

We summarize here a simplified theory of operation of the magnetometer (Ledbetter 
et al. 2008) to aid interpretation of our results. We use the simplifying assumption 
that the absorption linewidth is broad (18GHz/(atm buffer gas)) compared to the 87 Rb 
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Figure 1. (a) Laser beam path schematic. The 795 nm pump laser is linearly polarised 
by a polarising beam splitter (PBS) and then circularly polarised by the quarter 
waveplate (A/4). The 780 nm probe is detuned from the D2 resonance and linearly 
polarised. Both lasers are conveyed from their sources to the optics via polarisation 
maintaining (PM) optical fibres. The half waveplate (A/2) is used to balance the 
photodiode (PD) difference signal from the Wollaston polarising beam splitter (W). 
(b) A simplified diagram of the first three fine-structure 87 Rb energy levels. 



ground state hyperfine splitting (6.8 GHz), though in our cells the hyperfine splitting 
and pressure broadening are nearly equal and the model will be approximate. If the 
magnetometer is operated with sufficiently high Rb density and sufficiently low total 
magnetic field, spin-exchange collisions preserve the coherence of the transverse spin 
(Alfred et al. 2002) and three relaxation processes dominate the polarization relaxation 
rate T, 

r = r sd + r pr + r, (i) 

where R is the optical pumping rate (the rate that an unpolarised atom would absorb 
pump photons), r sd is the spin-destruction rate from collisions and r pr is the absorption 
rate of probe laser photons. 

Assuming the Rb atoms are in spin-temperature equilibrium (Anderson et al. 1960) 
and neglecting diffusion, the polarisation has a steady-state solution for static magnetic 
fields (Ledbetter et al. 2008) 

p _ p / + n a .r2 z \ ^ 

where ft = 7B is the Larmor frequency from a magnetic field B and 7 = g s fiB/h~- P 
is the Rb polarisation in the coordinate system of figure [T^a). The approximations on 
the right hand side of ([2]) hold for small magnetic fields such that |0| < T. At typical 
T « 2tt x 1000/s, this limits B to < 30 nT. 

2.2. Signal size and sensitivity 



The Rb polarisation modifies the susceptibility of the atomic vapor, which changes the 
index of refraction for a + and cx_ circularly polarised light. For linearly polarised light, 
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this causes a rotation of the polarisation angle. In our case, the probe beam is linearly 
polarised in the y — z plane and is rotated in this plane by (Ledbetter et al. 2008) 

1 



-r e fcnl- 



"0 



-P 



(3) 



where r e is the classical electron radius, / is the transition oscillator strength, c is the 
speed of light, n is the Rb atomic number density, I is the optical path length, v — u 
is the probe laser detuning from the D2 resonance, and Av is the pressure broadened 
absorption linewidth of the transition. For reference, using our experimental parameters 
(see section [6]), evaluates to <fi/B y ps 0.1/xRad/fT. For perspective, the differential 
photocurrent is AI = 2P e(f>, with probe laser power P and a photodiode power to 
current conversion e, is around 250pA/fT. 



3. Apparatus and design considerations 

3.1. Apparatus description 




Photodiodes 



Mirror 

Thermal 
Insulation 




(a) 

Figure 2. (a) Single magnetometer unit (left) with view inside the thermal insulation 
(right). Square tri-axial nulling coils are wrapped around the insulation, (b) (bottom 
insert) Four channel magnetometer array, with 7 cm channel spacing. This array is 
mounted on an existing SQUID gantry top-down, as in (b)(top). In human subject 
trials, the subject lays on a bed and the array is centered over the subject's heart by 
equalizing the signal in all four channels in real time. The x, y, and z axes in the 
text refer to magnetometer body axes. Here, all y are aligned and correspond to the 
sensitive direction for all magnetometers. 



Design choices in this atomic magnetometer are meant to minimize noise consistent 
with maximizing suitability for biomagnetic measurements. Figure [21(a) shows a single 
magnetometer unit. It consists of three fixed optics tubes with all the optics from the 
schematic in figure (Ha), a Pyrex 87 Rb vapor cell (Triad Technologies) press-fit inside 
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a boron-nitride heatsink thermally contacted to resistive-film heaters (Minco) and a 
thermistor, all surrounded by thermal insulation (Aerogel). A plastic scaffold clamps 
the heatsink, optics tube, and two custom-made plastic mirror mounts in place. Single- 
turn triaxial coils are wrapped around the thermal insulation to locally null the magnetic 
fields in the detection volume and provide calibration and control fields. Polarisation 
maintaining fibres connect the optics tube to the pump and probe lasers. The pump 
laser is a distributed feedback type with output power 80 mW (Eagleyard Photonics), 
detuned | Az/ pump | < 20 GHz from the Dl resonance. The probe laser is a 20 mW tunable 
external cavity diode laser (Sacher Lasertechnik Group), detuned ~0.2nm from the 
D2 resonance. The two lasers are each split equally into four separate fibres with an 
integrated polarisation maintaining splitter (Oz-Optics). The result is a set of individual 
magnetometer channels that are unconstrained with respect to one another, due to local 
fibre coupled optics and individual magnetic field control at each cell. 

The lasers are mounted on a portable breadboard which, along with all the coil 
control electronics and monitoring equipment, fit on a single cart. A PXI based FPGA 
data acquisition unit (National Instruments) is used to apply calibration fields and 
digitize magnetometer signals. The PXI crate also houses power supply cards to control 
the vapor cell heaters. This modular design ensures the entire apparatus is easily 
transported between the MSR and our testing lab in a different building. 

Material selection was aimed at minimizing Johnson noise caused by the use of 
electrical conductors. Each magnetometer unit shown in figure EJja) contains only a few 
conductive pieces, namely the heaters, thermistor leads, optical fibre terminations, and 
the photodiode leads. 

The vapor cell is 1 x 1 x 5 cm Pyrex rectangular container with isotopically pure 87 Rb 
and buffer gas composed of 50 Torr nitrogen to prevent radiation trapping (Happer 1972) 
and 760 Torr He at the time of manufacture. Subsequent measurement showed slow loss 
of buffer gas pressure indicating He loss consistent with diffusion through the Pyrex cell 
walls (Walters 1970). We estimate the data presented here was taken with a residual 
He pressure of ~200 Torr. The cells were operated at a temperature of 140-180 C. Each 
magnetometer requires 4W of heater power in thermal equilibrium, corresponding to 
electrical currents near 0.5 A, large enough to produce potentially disruptive magnetic 
fields and gradients. We manage these using matched heater pairs with oppositely 
aligned magnetic fields and a boron nitride heatsink to spatially separate the heaters 
from the detection volume of the vapor cell. Furthermore, the heater currents were 
modulated at 100 kHz so these magnetic fields and their gradient fields would average 
to zero on the time scale of MCG measurements. There is no detectable difference 
between the magnetic noise level with the heaters on or off. 

Figure EJ^b) shows a photograph of the array mounted to an existing gantry for a 
previously installed commercial SQUID system. Human subjects lie on a bed under the 
array. The array can be translated, rotated, and tilted using the gantry controls. For 
convenience in these initial tests, our measurements were made with planar square array 
and channel spacing of 7 cm. The minimum planar array spacing for all four elements is 
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4.5 cm. More magnetometers can be added, with 7 cm spacing in the vertical direction 
of figure [2](b) (bottom), while the optics tubes require the spacing in the horizontal 
direction to be >10cm. In our setup, the four elements operate independently and 
could be tilted or translated with respect to one another. It should be straightforward 
to implement non-planar geometries. 



3.2. Modularity 

DC fields in the MSR, 10-50 nT, decrease the fundamental sensitivity according to (J2j) 
and require nulling. Figure |3^a) shows a sweep of B y and a fit with R = 400 ~ T to 
(|2J) with Q X ,Q Z = 0. By must be nulled to better than \B y \ < 2nT to remain in the 
sensitive linear region of the response. 



2-i 




By (nT) By (nT) 

(a) (b) 

Figure 3. (a) Magnetometer signal as a function of applied field (dots) and a fit to 
an amplitude scaled version of ([2]l, (b) A simulation of P x for multiple magnetometers 
in the presence of the measured magnetic field gradients from the zero values at each 
channel. With only a single set of tri-axial coils, these are the minimum gradients in the 
shielded room, and represent the best-case operating point for the four magnetometers. 
In this case, a compromise value for B y must be chosen, reducing the usable operating 
range. This illustrates the need for local control over the magnetic fields at each cell. 



If the MSR magnetic fields were sufficiently uniform, we could use a single set of 
tri-axial coils with side length a and array channel baseline spacing d such that a 3> d 
to simultaneously cancel the magnetic fields at all magnetometers. Having tried this, 
we have found two problems with this approach that resulted in the choice of small 
individual cancellation coils for each magnetometer. 

First, residual MSR magnetic field gradients are too large on the scale of the spacing 
of our magnetometers to allow a uniform nulling field to be used. Figure [3](b) shows 
the response predicted by ([2]) for four magnetometers in the presence of measured MSR 
magnetic field gradients and R from the fit in figure [3]^a). B x , B y , and B z are the 
differences from the measured zero values at each magnetometer channel. As shown in 
figure EJ^b), these gradients cause a decrease in both the single channel sensitivity and 
the effective multichannel B y operating range (highlighted in yellow), which has been 
decreased by a factor of 4 to \B y \ < 0.5 nT. 
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Second, and most importantly, the pumping laser can cause a light shift (or AC 
Stark shift) due to the pump laser interacting differently with different Zeeman sublevels 
through the Rb vector polarisability (Savukov & Romalis 2005, Budker & Romalis 2007). 
This causes an effective magnetic field (Appelt et al. 1999) 

fl LS « RAs, (4) 

where A is the ratio of the laser detuning from atomic line centre to absorption linewidth, 
and s is the photon helicity, ±1 along the pump propagation direction for circularly 
polarised light. 

At high optical depths corresponding to high temperatures, the pump beam 
is absorbed before polarising the entire cell length unless it is detuned A ~ 1 from 
resonance, resulting in a nonzero J^ls- The magnitude | Ols I is roughly the same from 
cell to cell, but the direction with respect to the null coils depends on the orientation of 
the magnetometer and the direction of s. For our laser parameters, we expect the pump 
light to cause effective fields of magnitude ~50nT. Orienting all the magnetometers 
to have the same light shift with respect to the external coils mitigates the problem, 
though the nonuniform light shift within each cell due to the spatial intensity profile of 
the laser results in a broadening of the magnetometer dispersion (Appelt et al. 1999). 

The light shift presents a problem for the prospects of differing individual channel 
orientation and tilt, or nonplanar arrays. There are also practical problems associated 
with large coils, such as limited magnetometer adjustment and complications with 
calibration fields when the magnetometers are not aligned with the coil axis. These 
problems are all solved by using individual sets of coils for each magnetometer. 

4. Methods 

4-1. Residual Field Nulling 

We find the DC null fields for a single magnetometer as follows (Seltzer & Romalis 2004, 
Li 2006): 

• sweep By and set the operating point in the centre of the dispersive curve, where 
the magnetometer is most sensitive. 

• apply a 50 pT oscillating field B sm(u>t)x at 10-30 Hz and adjust B z until there is 
no response to B , ensuring B y remains at its most sensitive operating point. 

• apply B sm(uut)z and adjust B x to null the response, again ensuring B y is adjusted 
to maintain maximum sensitivity. 

• repeat until residual response to B x and B z is minimized. 

For the four- magnetometer array, we repeat this procedure for each individual unit. 
Iterating the procedure two times eliminates any offsets of the local field at the first 
magnetometer caused by adjusting the fields in the later ones to below the adjustment 
procedure sensitivity. 
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There is no cross-talk between our four channels because the magnetometers are 
fundamentally passive devices. Non-uniform light shifts from the Gaussian pump laser 
intensity profile cause non-uniform effective magnetic fields as in PJ . Due to the residual 
Q z term in fl2]), each magnetometer retains some sensitivity to B x . This is typically a 
correction of < 10% and its effect is ignored in what follows. 

As discussed in section 13.21 the magnetometers must be within approximately 1 nT 
of thier null point in order to maintain uncompromised sensitivity. Background field 
variations in the shielded room are generally smaller by an order of magnitude or more 
on the timescales of of a 30 s MCG measurement. Long-term variations can take the 
magnetometers outside of the optimal measurement range and are easily compensated 
between measurements. 



4-2. Magnetometer sensitivity and noise floor characterization 




Figure 4. Measured magnetic noise (solid black) of a single characteristic 
magnetometer channel on a magnetically quiet day in the MSR. (Inset) Amplitude 
response x(w) (open circles) and phase response 4>(uj) (dots) of a single magnetometer 
to a calibrated applied signal. The fits to x(w) and 0(w) are used to calibrate the 
all noise PSD and time-series measurements presented here. The horizontal line is at 
6 fT Hz -1 / 2 , the expected magnetic noise level based on table[2] The no-pump effective 
noise level (dashed blue) is 1-2 fT Hz^ 1 / 2 , consistent with our estimates for photon shot 
noise. 

Once the DC null fields have been determined, each magnetometer is simultaneously 
characterized by applying a set of 10-200 pT magnetic fields at twenty different 
frequencies ranging from 2-200 Hz, as in figure H]^ inset), providing a frequency dependent 
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calibration for each magnetometer signal S, such that S = x{^) el<l ^ w 'B y . We then take a 
short (~6 s) noise measurement detecting the residual response of the magnetometer to 
ambient fields and other noise sources. Each magnetometer is sampled at 20 kHz with 
a two-pole lowpass filter at 10 kHz, using a 16-bit digitizer. The power spectral density 
(PSD) of this signal is used to characterize the noise floor as 

5B y/pSD 

producing the measurement in figure HI We note that while the MSR floor is constant 
over time, sinusoidal residual magnetic fields often dominate the spectrum. Figure H] 
shows a measurement with a relatively quiet background, displaying only a few such 
peaks near at 6.7, 50.5, 60, 101, and 120 Hz. We suspect these peaks come from nearby 
air handling mechanics and powerline noise. Though this makes real time analysis of 
the magnetometer signals difficult, they tend not to interfere with MCG measurement 
because they are narrowband and easily identified. 

To obtain an estimate of the non-magnetic noise level of the magnetometers, we 
block the pump beam and measure the PSD once again. In this case, the device should 
not be sensitive to magnetic fields. This is non-magnetic noise that is transformed 
into an effective noise floor through the same calibration used in (jSJ). The residual 
noise is from photon shot noise, spin-projection noise, and non-magnetic technical noise. 
Figure m inset) shows an example of the single magnetometer calibration parameters 
x(o;) and <j)(uj). The blue dashed curve of Figure H] shows the non-magnetic noise level. 
All magnetometers operated with a baseline sensitivity < 11 fT Hz -1 / 2 and a probe-noise 
limited sensitivity at least a factor of two below the measured magnetic field sensitivity. 

5. MCG 

We obtained MCG from 13 adult human subjects. The experimental protocol was 
approved by the institutional human subjects committee prior to commencement of the 
study and informed consent was obtained from all subjects. 

Each measurement consists of a 30 s MCG recording from each magnetometer 
sampled with the same parameters described above. For data analysis, the magnetic 
field in the time-domain is obtained by a deconvolution of the magnetometer response 
(measured in the calibration) from the measured signal. The signal is then down- 
sampled to 1024 Hz. Low frequency drift is removed using a wavelet detrending 
algorithm. 

Figure E^a) shows a processed segment of 30 s MCG recordings from the four 
channels. Figure EJ^a) shows a reference measurement taken just after the data presented 
in figure O but without a subject. This measurement was taken on a different day than 
figure|4]and there are clearly more peaks, mostly at 4.55 Hz and its harmonics, along with 
a prominent 6.7 Hz peak and powerline harmonics. To obtain the traces in figure 151(a), 
the magnitude and phase of these non-biomagnetic peaks is determined from the 
reference, and then a matching sinusoid is subtracted from the signal channels, resulting 
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Figure 5. (a) Four-channel adult MCG with large signal interferences subtracted, (b) 
Shows the same set of MCG after averaging each channel. 
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Figure 6. (a) MSR noise background for the MCG data. The peaks are actual 
magnetic artifacts, likely from building air-handling units, and comparison with 
figure 21 obtained on a different day, illustrates the day-to-day variation in the MSR 
background fields, (b) MCG signal (black) and MSR noise background (red dots) after 
the signal processing described in the text. The non-biomagnetic peaks in (a) are easily 
removed. 
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Table 1. QRS peak amplitudes of this work. QRS dipole moments are estimated 
to compare with previous works. The distance range quoted here is for the range of 
individual magnetometer distances to the heart. 





(Nousiainen et al. 1994) 


(Knappe et al. 2010) 


This work 


r (cm) 


8-10.5 


5 


5.5-10.8 


•Bqrs (pT) 


(not reported) 


100-150 


40 ± 16 


m (/xAm 2 ) 


0.76 ±0.36 


0.4-0.6 a 


0.60±0.48 a 



Estimated from -Bqrs and r. 



in the PSDs shown in figure [6(b). Figure [5(b) shows the averaged traces. The P, QRS, 
and T components are well resolved. We emphasize these are magnetometer signals 
(rather than gradiometer signals), and contain any unfiltered background magnetic noise 
in the MSR as well as the MCG signal. 

We estimate the minimum distance to the skin for each magnetometer of around 
1 cm. Knappe et al. (2010) reported their sensor head volume was located 5 mm from 
the skin or 5 cm from the heart, whereas the SQUIDS used for comparison were 7.5 cm 
from the heart. With our channel spacing (7cm), we estimate each magnetometer is 
7.6 cm from the heart when the array is centered above it. 

A statistical summary of the QRS peaks of 13 subjects is presented in table HJ 
along with an estimation of the QRS dipole amplitude m. We use the approximation 
-Bqrs ~ m/r 2 - 6 (Nousiainen et al. 1994) to account for distributed source effects. To 
calculate m, we require the distance from the heart to each magnetometer. We minimize 
the standard deviation of m across the four magnetometers as a function of the heart 
position with respect to the center of the magnetometer array, which can then be used 
to calculate the dipole amplitude for each patient. The result from this work in tabled] 
is a weighted mean of the 13 measured subjects. 

In addition to adult MCG measurements, we also simulated MCG using a head 
phantom (Biomagnetic Technologies), which is a spherical shell filled with saline solution 
and five different current dipoles that can be driven individually. The head was placed as 
close as possible to the array centre and aligned so that the detected signal from driving 
a dipole was maximized. We drove the phantom with a test waveform (Goldberger 
et al. 2000) shown in figure [7(b). 

We used the calibrated magnetometers to adjust the QRS peak height driven by 
an amplitude scaled version of the test waveform. Figure [7(a) shows 28-beat averages 
for all four magnetometer channels, after processing the signals in the same way as the 
MCGs. The detected QRS amplitude is about 2pT, and the averages for each channel 
faithfully represent the field expected from the driving waveform. 
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Figure 7. (a) 28-beat average response from the phantom driven with the waveform 
shown in (b), with magnetic field offsets applied to aid visualization. Before averaging, 
the signals are processed in the same way as the MCG above, by removing only narrow 
frequency components. 

Table 2. Magnetometer Noise Budget. Two results are quoted for the photon shot 
noise. The basic model is (0 where P z , T pr , and T are all averages over the laser 
interaction volume. In the detailed model, these quantities are calculated at each 
position within the cell and then integrated to recover the photon shot noise floor. 



Noise Source 


Level (fTHz -1 / 2 ) 


MSR Johnson noise 


5 


Electrical current noise in nulling coils 


3.9 


87 Rb Johnson noise 


0.9 


Photon shot noise (basic model) 


0.3 


Photon shot noise (detailed model) 


2 


Spin-projection noise 


0.03 


Total 


6.2 



6. Noise sources and characterization 

There are several different sources of unwanted signal in our magnetometer. These fall 
into three general categories. One is technical noise, including electrical pickup noise or 
probe movement on the detector photodiodes. The second is magnetic signal other than 
the desired biomagnetic signals, and the third noise category is fundamental quantum 
noise associated with the magnetic measurement. The last two terms are the dominant 
noise sources and we estimate the noise floor in table EJ 
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6.1. Magnetic noise 

Magnetic noise generally includes any magnetic source that falls within the frequency 
band of interest and limits the achievable SNR of the desired signal. Clear examples 
can be seen in figure E^a). 

The single channel magnetometer noise floor in figure BJb) and figure [6] is limited by 
Johnson noise in nearby electrical conductors. These conductors carry thermal currents 
which generate white-noise profile magnetic fields. The most important noise sources 
in our case are likely the mu-metal walls of the MSR and 87 Rb thin-film remaining on 
the Pyrex cell walls. We estimate these contributions from the calculations of Lee & 
Romalis (2008). The 2m cube shielded room has a noise floor of around 5fTHz _1 / 2 . 
We roughly model the 87 Rb film as a cylinder with a 1 cm radius and 10 nm thickness, 
about 5 mm from the detection volume, which would have a 0.9fTHz _1//2 noise level. 

We also measured the electrical current noise in the circuitry used to drive our 
nulling coils using a transimpedance amplifier. This contributes 3.9fTHz -1 / 2 to the 
magnetic noise level. 



6.2. Photon shot noise 

Though most SERF magnetometers are limited by shield Johnson noise, there are two 
fundamental quantum noise sources which limit the maximum theoretical magnetic field 
resolution: spin-projection noise and photon shot noise. We briefly discuss the results 
in Ledbetter et al. (2008) as they pertain here to photon shot noise. We estimate 
spin-projection noise to be a negligible contribution to the overall noise budget. 

Photon shot noise occurs because the polarisation of the probe laser is determined 
by counting the difference in the number of photons arriving at the two different 
photodiodes in figure QJa). The resulting uncertainty in the magnetic field is 



y P z ^Y pr OD Q nVt' 

where OD = 2r e fcnl/ Aj is the probe optical depth on resonance. 

The parameters used in our calculations were R = 2575/s, r s d = 60/s, r pr = 155/s, 

n = 7.8 x 10 13 /cm 3 , pump waist w 0:PU = 2.6 mm, and probe waist w 0:PT = 1.1mm, and 

interaction volume V = 7ri>q 2iOo,p U - Equation (jBJ) estimates a photon shot noise level 
of O.SfTHz- 1 ^ 

We have found using this calculation with spatially averaged laser parameters 
estimates a lower photon shot noise limited sensitivity than we actually achieve in our 
experiments. A more detailed (Walker & Happer 1997) model approximating the pump 
and probe laser propagation through the cell and their Gaussian spatial profiles, gives 
a more realistic 2fTHz -1 / 2 . 

For comparison, an experimental measure of the photon shot noise can be obtained 
using the measured differential photocurrent per unit field AI p /AB and an effective 
photocurrent shot noise where I p is the photocurrent from a single photodiode. The 
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photon quantization noise is 5I p /y/Hz = ^/2eJ p . The effective magnetic field noise is 
then 

SB = 5I p /VTti 

AI P /AB' U 

For a typical measured AI p /AB = 8 /iA/nT operating with all four channels 
and a total probe power of 1.67 mW, SB/Vm, = l^fTHz" 1 / 2 from in reasonable 
agreement with the more detailed calculation. 

7. Conclusions and future work 

We have presented a modular four-channel atomic magnetometer array operating in 
the SERF regime. The adjustable 7 cm channel spacing used for these measurements 
is much larger than previously presented biomagnetic measurements using SERF 
magnetometers. The single channel sensitivity ranged from 6-11 fT Hz -1 / 2 , dominated 
by magnetic noise from the MSR. The modularity of the magnetometers, in principle, 
allows flexible channel positioning in an array, despite large light shifts in different 
directions in the room reference frame. 

Our sensitivity is limited by real magnetic noise. In this case, parametric 
modulation (Li et al. 2006) can allow simultaneous measurement of B x and B y , at 
the cost of a small (1/2-1/4) decrease in signal size. Future work will aim at using 
parametric modulation to measure two components of the magnetic field with each 
modular unit. 

Detection of fMCG, an application which requires high sensitivity but a modest 
number of channels, should be possible with our array. However, the number of channels 
is insufficient to allow the effective use of spatial filters, such as beamformers and 
independent component analysis, which are widely used to remove maternal interference 
and to improve the signal-to-noise ratio of fMCG recordings. With a modest reduction 
in size, the design can be extended to higher channel counts. This will be essential for 
constructing dense arrays suitable for MEG and adult MCG. 
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